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Abstract 

This paper focuses on some simple models of limit order book dynam- 
ics which simulate market trading mechanisms. We start with a discrete 
time/space Markov process and then perform a re-scaling procedure lead- 
ing to a deterministic dynamical system controlled by non-linear ODEs. 
This allows us to introduce approximants for the equilibrium distribution 
of the process represented by fixed points of deterministic dynamics. 

1. Introduction. The underlying Markov process 

In what follows, LOB stands for the limit order book, a trading mechanism 
adopted in many modern financial markets. For a detailed description of some 
common LOB models and their applications, see [T] and references therein. (Al- 
though our models differ in a number of aspects.) One of challenging problems 
is to determine factors attracting or repelling various market participants. 

This paper explores a new approach to the analysis of LOB dynamics where 
the parameters of the original random (Markov) process are re-scaled, and a 
limiting dynamical system emerges, with a deterministic behavior described by 
a system of non-linear (ordinary) differential equations. A similar approach is 
commonly used in the literature on stochastic communication networks; see, 
e.g., the paper [2] and its sequels (in particular, [3]). In the current paper we 
consider a simplified model where a number of technically involved issues are 
absent. We also omit proofs, referring the reader to forthcoming publications 
beginning with [4]. 

The rationale for the models below is as follows. We consider a single- 
commodity market where prices may be at one of N distinct levels (say, ci < 
C2 < ... < Cat, although the exact meaning of these values is of no importance). 
The market is operating in discrete times 0, S, 2(5, .... At a given time tS, 
t = 0, 1, . . ., there are bi{t) traders prepared to buy a unit of the commodity at 
price Ci and Si{t) traders prepared to sell it at this price, which leads to vectors 

bit) - {hit),...,bM{t)), s{t) = {si{t),...,SN{t)) eZ^. (1) 
Here Z+ = {0, 1, . . .} stands for a non-negative integer half-lattice and for 
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the non-ncgativc integer A^-dimensional lattice orthant. The pair {b{t), s{t)) 
represent a state of a Markov process U{t) that will be the subject of our 
analysis. 

If bi{t) > Si{t) > then each of the sellers gets a trade with probability G 
(0, 1) and leaves the market, together with his buyer companion. Therefore, 
both values bi{t) and Si{t) decrease by a random number n = 0,1, ... ,Si{t) 
with the binomial probability. A seller among Si — n who did not get the trade 
either (i) quits the market with probability pq € (0, 1) or (ii) moves to the price 
level Ci-i with probability pu € (0, 1) or (iv) remains at the same level with 
probability 1 — pq —pu- (One can think that for this seller a random experiment 
is performed, with three outcomes.) Similarly, a buyer among bi(t) — n who did 
not get the trade quits the market with probability pg € (0, 1) or moves to the 
price level c^+i with probability pu G (0, 1) or remains at the same level with 
probability 1 — pQ — pm- (Assuming that pQ + pu < !•) 

Symmetrically, if Si{t) > bi{t) > then each of the buyers gets a trade 
with probability px £ (0, 1) and leaves the market, together with his seller 
companion. The remaining traders at the price level Cj proceed as above. 

In addition, at every time point tS a random number of new buyers arrive 
and position themselves at the price level ci; it is distributed according to a 
Poisson law with mean Aj, > 0. Similarly, at every time td a Poisson random 
number of sellers arrive and take a position at price level cjv; the mean value of 
this variable is Ag > 0. 

All described events occur at each level independently. This generates the 
aforementioned Markov process {C/(f)} with trajectories {{b{t), s{t))}, t € 1^+. 

Theorem 1. V values of parameters A^/s, Pq/m andpx, the process {U{t)} 
is irreducible, aperiodic and positive recurrent. Therefore, it has a unique set 

of equilibrium probabilities n = and V initial state U{0) 

(deterministic or random), the distribution of the random state U{t) at time t 
converges weakly to n as t oo: 

lim P(c/(t) = (6,s)) =7r(6,s). (2) 



2. Scaling limit 

The explicit form of the equilibrium distribution tt of process |C/(t)} (and 
even probabilities of transitions (6, s) (6', s')) are too cumbersome to work 
with. This fact makes it desirable to develop efficient methods of approximation. 
In this paper we focus on one such method based on scaling the parameters of 
the process (including states and time-steps). 

The re-scaling procedure is as follows: we set 

7 "Q Q!M . Ab . As 

Pt = y, PQ = Pm = -j-, Ab = — , As = — , (3) 
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where 7 > 0, ag > 0, au > 0, Ab > and As > are fixed and L — > 00. In 
addition, we re-scale the states and the time: pictorially, 

6i £j M 

^'^L' y^^L' 

Formally, denoting the Markov process generated for a given L by C/*-^-', we 
consider the continuous-time process 

V^^\t) = \u^^'^{\tL/5'\), T>Q, (4) 

where \a\ stands for the integer part of a > 0. 

Set: R+ — (0,00) (a positive half-line), then M.^ is a positive orthant in N 
dimensions. Suppose we are given a pair of vectors (a;(0), y(0)) £ M;'!' xM;!^ where 
a;(0) = (xi(0),...,XAr(0)), y(0) = (yi(0), . . . , yAr(O)). Consider the following 
system of first-order ODEs for functions Xi = Xi{T) and j/^ — yi{T) where r > 
and 1 < i < TV : 

ii = Ab - (^aq -I- aM^a;! - 7min 

Xi ^ aMXi-i - [uq + aMjXi - ^ uAn [x^, y^] , 1 < i < iV, 
m = OLMVi+i - f aq + olm) Vi - 1 min [a;i, , \ <i < N, 
VN ^ - (aq + aujvN - ^rain [xn.Un], 

with the initial date 2:^(0), yi{0), 1 < i < N . The fixed point (a;*,y*) of system 
^ has X* — (x*, . . . , x^) and y* = {yl, . . . , y^) where x* and y* give a solution 
to 

Ab = (aq + aMjxl+"f min [x* , yl] , 

aM<_i = (aq + aM)a;* +7min [x*,y*], 1< i < N, 
} \ 

aMy*+i = [o'Q + o^Mjyi + 7 min [x*,y*], 1 < i < N, 
As = (aq + aM^y^r + 7min [x^,y;^]. 

Both systems ^ and © are non-linear. However, the non-linearity 'disappears' 
at a local level which greatly simplifies the analysis of these systems. 

In Theorems 2 and 3 below, we use the distance generated by the Euclidean 
norm in x R^. 

Theorem 2. V initial date (a;(0),y(0)) £ R^ x R^ there exists a unique 

solution ^{x{T),y{T)), t > o| to system (5). For this solution, {x{T),y{T)) G 

R^ x R^ V r > 0. As T 00, the solution approaches a fixed point, which 
yields a unique solution to system (6): 

dist[(a;(r),y(T)),(a;*,y*)l ^0. (7) 
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Theorem 3. Suppose that the re-scaled initml states —U{0) tend to vector 

L 



(a;(0), y(0)) G x in probability: V e > 0, 



lim P dist 



-?7(0),(a:(0),y(0)) 



> e = 0. 



(8) 



Then, V r > 0, the process |v^^-'(t), t G [0,T]| converges in probabiUty to 
the solution |(a;(r), ^(t)), < t < r|. That is, V e > 0, 

^lim P (sup |dist V'-^\t), {x{T),y{T)) : < r < t| > = 0. (9) 
III particular, if x{0) = x* and y(0) — y* then 



< r < T 



0. 



(10) 



lim P (sup (dist \v^^\T),(x*,y' 

Moreover, if process {U{t), t G Z+} is in equilibrium then Eqn (10) holds true. 

3. Fixed points in the scaling limit. Concluding remarks 

The approximation developed in Theorem 3 calls for an analysis of solutions 
to (6). As follows from the middle equations in (6), 



Lemma 4. The fixed-point entries satisfy 

xl> .. . > x*ff and yl < . . . < y%. 



(11) 



Consequently, the parameter space formed by 7, ag/M; ^^d. Ab/s parti- 
tioned into open domains where one of the following generic patterns persists: 
(i) x% > ylf, (ii) X* < 2/*, and (iii) x* > y* for i = I, . . . , £ and x* < y* for 
i = £ + 1, . . . , N where I < £ < N . In each of these domains system ([6]) is linear. 



Lemma 4 allows us to develop simple algorithms for calculating the fixed 
point (a;*, y*) and analyze the character of convergence in ([7]). 

A particular algorithm for calculating {x*,y*^ is based on the following 

recursion. Set xf^ = and let yf^ be the solution to the third and the forth 

equations ^ with y^^ = As/(aQ + ckm)- Next, let (a^C"), yC^'), fc = 1,2,... be 
the solution to the system 



Ab = ("Q + ctM^ xj''-' + 7 min [xf^ , yf ^''] , 
ayixiJi = ( aq + aiAjX] ' + 7 mm [xl ' ,yl '\, I <i< N, 



amy, 

As = (aq + auj + 7 min [x^ , y^] 
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Lemma 5. The inequalities x^^"^ > xf^ yf^ < yf' hold true 

(k) . . (k) 

y i,k > 1 and, values xl are uniformly bounded. Therefore, 3 lim , and 
lim and these limits satisfy the system 

We conclude with the following remarks. 

1. The current set-up admits straightforward generalisations to the case 
where parameters 7 and ckq/m depend on z, < i < A'^ and on the trader 
type (b/s). A more complicated case emerges if parameters Ab/s become state- 
dependent. 

2. There are several forms of convergence for which the assertion in Theorem 
3 holds true. The dynamical system ^ itself gives rise to a limiting process 
with interesting properties. 

3. Another valid approximation for process U{t) is a diffusion approximation 
working on a different scale from that in (jSj. 

These topics are subject to forthcoming research. See [i\ and subsequent 
publications. 
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